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We analyze the semiclassical Einstein equations for quantum scalar fields satisfying modified 
dispersion relations. We first discuss in detail the renormalization procedure based on adiabatic 
subtraction and dimensional regularization. We show that, contrary to what expected from power 
counting arguments, in 3 + f dimensions the subtraction involves up to the fourth adiabatic order 
even for dispersion relations containing higher powers of the momentum. Then we analyze the 
dependence of the trace of the renormalized energy momentum tensor with the scale of new physics, 
and we recover the usual trace anomaly in the appropriate limit. We also find selfconsistent de 
Sitter solutions for dispersion relations that contain up to the fourth power of the momentum. 
Using this particular example, we also discuss the possibility that the modified dispersion relation 
can be mimicked at lower energies by an effective initial state in a theory with the usual dispersion 
relation. 
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I. INTRODUCTION 

The expansion of the Universe can act as a cosmological microscope. Scales of interest today could have 
been sub-Planckian at the first stages of inflation, and it has been therefore speculated that there could be 
signatures of trans-Planckian physics in the evolution of the universe and/or in the inhomogeneities of the 
CMBR [H, Q . It is of course very difficult to address these issues without knowing the physics at Planck 
scale. A plausible approach is to test the robustness of inflationary predictions under different changes that 
could be attributable to the unknown new physics. For example, it has been argued that one of the effects 
could be the modification of the dispersion relations for the modes of the quantum fields with momenta 
larger than a given scale M c , as suggested by loop quantum gravity [|[ or by the unavoidable interaction 
with gravitons [1] . It is therefore of interest to study the consequences induced by such modifications into 
the dynamics of the scale factor of the Universe. This is the main purpose of the present work. We will 
consider a quantum scalar field with a modified dispersion relation on a Robertson Walker background, and 
analyze the Semiclassical Einstein Equations (SEE) in order to test whether the physics at very high scales 
may leave an imprint at lower scales or not. 

The source of the SEE is the mean value of the energy momentum tensor of the quantum fields, {T^J), 
which is formally a divergent quantity. For quantum fields with usual dispersion relations, the covariant 
renormalization procedure is very well known [5|, 0, 0] ■ For Robertson Walker metrics in n dimensions, one 
can compute the energy momentum tensor in the so called adiabatic approximation [j| [l^ , and define the 
renormalized one as 

(T^)ren = (T^) - (T^ ... - (T^ 2 ^ (1) 

where [n/2] is the integer part of n/2. (T^)^' denotes the terms of adiabatic order j of (T M „) (i.e., the 
terms containing j derivatives of the metric) and once regularized is proportional to a geometric tensor. 
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This procedure is equivalent to a redefinition of the gravitational bare constants of the theory. In order to 
give sense to the divergent quantities appearing in the equation above it is necessary to use a regularization 
method, a very useful one being dimensional regularization. In 3 + 1 dimensions divergences arise up to 
the fourth adiabatic order, and thus it is necessary to include in the gravitational part of the theory terms 
quadratic in the curvature. In toy 1 + 1-dimcnsional models, the fourth adiabatic order is finite, and the 
rcnormalization only involves the subtraction of the zcroth and second adiabatic orders. 

The SEE have also been considered for quantum fields with modified dispersion relations (Tll.[l2l|. However, 
these previous works do not include a proper treatment of the divergences, and the infinities are removed by 
simply neglecting the zero point energy of each Fourier mode of the quantum fields. This approach is not 
equivalent to a redefinition of the bare constants of the theory, and therefore is not fully justified. 

The extension of the adiabatic renormalization for generalized dispersion relations was first considered by 
us in Ref . [Til ] , where we computed the mean value of the energy momentum tensor up to the second adiabatic 
order. Power counting arguments suggest that, for dispersion relations containing fc 4 or higher powers of 
the momentum, it would be enough to subtract up to the second adiabatic order in 3 + 1 dimensions and 
up to the zeroth adiabatic order in 1 + 1 dimensions. However, a more careful analysis in 1 + 1 dimensions 
showed that this is not the case indeed, the quantities to be subtracted must be expressed in terms 
of geometric tensors in n dimensions, and the renormalization should be performed before the limit n — > 2 
is taken. In 1 + 1 dimensions G^ v is proportional to n — 2, since it results from the variation of the would 
be Gauss-Bonnet invariant in n — 2. Thus, the mean value of the energy momentum tensor is written as 
(Tfj, v )^ = c n G^ v where the constant c„ contains a pole at n = 2. Therefore, the second adiabatic order 
should also be subtracted in 1 + 1, although when explicitly computed it is a finite quantity. In this way 
one obtains the correct trace anomaly and the results are continuous in the limit M c — ► oo [Til ]. In this 
paper we will show that the same situation arises in 3 + 1 dimensions with the fourth adiabatic order. The 
tensor that results from the variation of the would be Gauss-Bonnet invariant at n = 4 is proportional to 
n — 4, and therefore (T^)^, although finite, contains a pole when expressed in terms of geometric tensors 
in n dimensions. The conclusion will be that the subtraction must also include the fourth adiabatic order, 
whatever the dispersion relation. 

As the calculations of the fourth adiabatic order are technically rather involved, for the benefit of the reader 
we would like to describe here with some detail the organization of the paper, our main conclusions, and the 
relation with previous works on the subject. In Section 2 we present the Lagrangian and energy momentum 
tensor of the scalar field with modified dispersion relations, as well as expressions for the geometric tensors 
in dimensional Robertson Walker spacetimes that will be used along the rest of the paper. Section 3 
describes the WKB mode functions for the scalar field, up to the fourth adiabatic order, and for generic 
dispersion relations. For the usual dispersion relation, the fourth adiabatic WKB mode functions have been 
computed previously in Ref. 1101] . while for generic dispersion relations, they have been computed up to the 
second adiabatic order in Ref. [l3[ . In Section 4 we use the WKB mode functions to construct the regularized 
{Tfw) i with j = 0, 2, 4 in n dimensions. These results generalize our previous ones for j = 0,2, and allow us 
to discuss one of the main points of this paper, which is the necessity to subtract up to the fourth adiabatic 
order in n = 4. The calculations arc relatively straightforward but require a lot of algebra, so we relegate 
the details to the Appendix A. In Section 5, and as a warm up, we analyze the renormalization of the stress 
tensor in n — 2. We compute explicitly the trace of the energy-momentum tensor in de Sitter spacctimc, 
and show that it reproduces the usual value as M c — > oo. Section 6 deals with the renormalization of (T^) 
in n — 4. Once more, as an example, the trace of the energy-momentum tensor is computed in de Sitter 
space for a massless conformally coupled quantum field, and the usual trace anomaly is recovered in the 
limit M c — ► oo. 

In Section 7 we present a concrete application of the formalism developed in previous sections. We first 
show that in de Sitter space, for any dispersion relation, there exists a one complex family of quantum states 
for which the unrenormalized energy momentum tensor is proportional to g^ v . One member of the family 
corresponds to a rcnormalizablc state, and for this particular state it is possible to obtain sclfconsistcnt 
de Sitter solutions, as explicitly shown for massless conformally coupled fields. We find that, for values of 
the cosmological constant smaller than a critical value, the SEE admit two different de Sitter solutions, a 
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perturbative one, close to the classical solution, and a nonperturbative one, with a very high curvature. 

In several previous works (l5l . [H| , it has been argued that the trans-Planckian physics could be taken into 
account by considering an effective field theory with usual dispersion relations, in which the effects of the 
new physics is encoded in the state of the field modes when they leave the sub-Planckian regime. In Section 
8 we discuss the possibility of simulate the modified dispersion relations by an effective initial quantum state 
at the level of the SEE, i.e. to discuss the backreaction effects. 

Throughout the paper we set c = 1 and adopt the sign convention denoted (+++) by Misner, Thorne, 
and Wheeler 



II. THE MODEL 



We consider a free quantum scalar field <f> with modified dispersion relation propagating in a curved 
space-time with a classical spatially flat FRW metric given by 

ds 2 = g^dx^dx" = -(u fl dx tl ) 2 + ± fll , dx^dx" = C(r))[-drj 2 + <%cfa;W] (2) 

where fi, v = 0, l...n — I (with n the space-time dimension), C x l 2 ir\) is the scale factor given as a function 
of the conformal time 77, the vector field = C 1 / 2 ^)^, and J_^= g^v + u^u v coincides with the spatial 
metric as defined by an observer comoving with u^. 

The classical action for the scalar field can be written as [ll| : 

S<jj = J d n x^J~ r g{C t j > + Ccor + C u ), (3) 

where g = det^g^), is the standard Lagrangian of a free scalar field 

Cf, = - i [g^d^d^ + (m 2 + £i?)</> 2 ] , (4) 
with R the Ricci scalar, C cor is the corrective lagrangian that gives rise to a generalized dispersion relation 

C cor = ~ M^VX^), (5) 

s,p<n 

with T) 2 § =J-n Va -Lit V 7 </> (where V M is the covariant derivative corresponding to the metric g^ v and -L^= 
g Xv -L^i/), and C u describes the dynamics of the additional degree of freedom u M whose explicit expression 
is not necessary for our present purposes. 

The generalized dispersion relation takes the form 



^ = fc 2 +C(7?) 



(6) 



where b sp are arbitrary coefficients, with p < s. 

The Fourier modes \k corresponding to the scaled field x — C'™ -2 "'/ 4 ^ satisfy 

with the usual normalization condition 

XkXk* ~ x'kXk = i- (8) 

Here primes stand for derivatives with respect to the conformal time 77, and in the conformal coupling case 
we have £ = = (n — 2) /{An — 4), while £ = corresponds to minimal coupling. 
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On the other hand, the gravitational action is given by 



So = 



16ttG 



(9) 



where i?^„p CT is the curvature tensor, R plJ = R p „ pll , and A, Gjv, a, P, and 7 are bare parameters which are 
to be appropriately chosen to cancel the corresponding divergences in the (T^) derived from S^. It is well 
known that in n = 4 dimensions, the quadratic terms in Eq.@ are necessary for renormalizing the effective 
theory in the case of scalar fields with the usual dispersion relation. Moreover, as we will show, they are also 
necessary for any of the generalized dispersion relations of the type given by Eq. Q . 
In the frame defined by the vector field w M = C 1 / 2 (ry)^, we can write the SEE as 



8nG 



(10) 



where (T^J) is the expectation value of the energy momentum tensor of the scalar field which in this frame 
satisfies {T^)- p = 0, G M „ = i? M „ — g^R/2 is the Einstein tensor, and 



H 



(2) = 



H pi/ — ri 9(ii/R>p8(7'yR> p ^ 2Rp u p^ f(J R v p A^R^y -\- "IR-^u -\- AiR^fjR v AR P Rpp^av- 



(lla) 
(lib) 
(11c) 



The particular combination 



Ha„ + H® - W W = H$ (n - 4), 



(12) 



which comes from the variation of the would be Gauss-Bonnet topological invariant at n = 4, will be very 
important for our discussion about the renormalizability of the theory in four dimensions. 

For the conformally flat metric that we are considering, the covariantly conserved tensors in Eq. (|ll[) are 
not independent, 



2fT$ + (n - l)[(n - 2)H^ - 4H$] = 0. 



(13) 



Therefore, it is enough to work with H^v and H^J , whose non trivial components are given by 



H 



C 

(i) _ 2(n-l) 



H (D = _ 



(n-10)(n-2) w4 



32 



C 



-H'H 2 



£(„-4)(n-6)-2 [± 



16 



= - 



(n - l)(n-2)(n - 3) 
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32C 
(n-2)(n-3) 
4C 



H 4 , 



, Q- 5) ^ 4 



-2 + ("- 4 ) 



("~4) 
64 



(n 2 - 13n + 28) 



(14a) 
(14b) 

(14c) 
(14d) 



with 



(1,3) 
11 



1,3) 
22 



if, 



(1,3) 

(n-l)(n-l) 



iy and W ee C"/C. 
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The nontrivial components of the Einstein tensor are 

Gn = G-22 = ■■■ = G(n-l)(n-l) = 

and the Ricci scalar takes the form 



(n-1) (n-2)„ 2 



4 2 
(n-2) 



R = 



(n-1) 
C 



The expectation value of the energy momentum tensor of the scalar field can be written as [2C 



+ £ 



c 

' (2 7 r v / C)(™- 1 ) 
(n-1) 



Vc("- 2 )/ 4 / 

/2 



2 2 



C / / * / *\ 2) i i2 

-(jiXkXk + XkXk ) - -Q2 2 — 



(T n ) = 



d^k^ I (1 _ \ c(n _ 2)/2 / 
(2^VC)("- 1 ) I \2 J V 



XA- 



fc 2 \ 



C(n-2)/4 



+ fGu|Xfc| 



1 / dfc 2 2 



IXfcl 2 -^(XfcXfc + XfcxT) +£^;(XfeXfc + Xkx'k) 



+ £ 



(n-2) f C" (8 - n) C' 2 



C 



4 C* 2 



to I 



(15a) 
(15b) 

(16) 



(17) 



(18) 



with Tn = T22 = ... = T( n _i)(„_i). Here n is the physical space-time dimension, and fi is an arbitrary 
parameter with mass dimension introduced to ensure that x has the correct dimensionality. 



III. THE WKB EXPANSION 



In order to compute the WKB expansion, we begin by expressing the field mode Xk in the well-known 
form 



1 



Xk 



: CXp 



W k (fj)dfj 



Vim 

Then, the expectation values of the stress tensor can be written as 



(19) 



(T v 



^TI-l 



C f dkk n - 2 ti n -" J [(W 2 )'] 



(27TVG) 



32W| 



.Gjrrn ^ (n - 1) 
? W k ? 2 



Wk 
2 

2V 



(n-2) 



2W k 



C'\n-2) , C"(W 2 ) 
16W fc C 2 



8CW| 



^ (n-2) , G'(1T 2 ) 
C 2 21T fe 



C 2FT fe 3 



(20) 
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11/ 




"1 , (n-2) 



2VF fe 
3 [(^ fc 2 )f (n 



C" 2 (n-2) C'(W 2 )' 



16IF fc C 2 



8CW| 



i)C(w?y 



cwl 



(21) 



where we have defined the factor fl^ = 2n ( - n - 1 ^ 2 /T[{n - l)/2] coming from the angular integration. 
Substitution of Eq.lO into Eq.® yields 



El 



3 



/2 N 



2 W| 



(22) 



where 0,1 = 



(£ — £„) Ci£. This non-linear differential equation can be solved iteratively for Wk by 



assuming that it is a slowly varying function of r\. In this adiabatic or WKB approximation, the adiabatic 
order of a term is given by the number of time derivatives. 

In what follows we will denote by ^W 2 the terms in W% of adiabatic order j. We obtain straightforwardly 
( 2 )W| replacing Wk by on the right hand side of Eq. (f2"2")) . 



(2)^2 



C" 

(€-£»)(n-l)l-^ 



191 

4 C 
16 C 2 



1_ ZJf dfc 2 
_ k 2 dwj 

1 ~~ UJ 2 dfc2 



(n - 6) C' 2 \ 
4 C 2 J 

lC" 2 fc 4 ri 2 a; 2 
4 C 2 w 2 d (fc 2 ) 2 



^[/ 2 - 4/ + 4(£ - £ n )(n - l)(n - 2)] + ^[/ 



where we have used the fact that uil = ufl/C depends on time only through the variable x 
the right hand side in terms of the function 



4(£-£ n )(n-l)], (23) 
= k 2 /C to rewrite 



ciln ui\ 
din x 



1. 



(24) 



and a dot means a derivative with respect to lnir. 

For n = 4 dimensions we will also need the fourth adiabatic order, therefore, we perform one more iteration 
in Eq. (|22|) . Discarding the higher order terms we arrive at 



^W 2 



where we have used that (uj' 2 .)' /to k 



-j^ {^W 2 [H 2 (2f if 2 ) 2H'f] 5^W 2 'Hf 2^W 2 "} , 



(25) 



2v /, ,2 _ _ Hf and ( W a)'7 w a = W 2 (/2 + f) _ n , f _ 



In what follows it will be relevant to know the dependence with k of the different adiabatic orders. From 
Eqs. (|22|) . |23|) and ([25]) . with the use of an inductive argument, it can be shown that the 2j— adiabatic 
order scales as u>t~ 2 -' . 
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IV. REGULARIZED ADIABATIC STRESS TENSOR 



In this Section we present the regularized adiabatic stress tensor, up to the fourth order, in n dimensions. 
Our goal will be to show that the different adiabatic orders are proportional to the conserved tensors g^ v , 

and Hpl that appear into the SEE. Therefore its subtraction will be equivalent to a redefinition of the bare 
gravitational constants of the theory. 

In order to find the tensorial structure of the different adiabatic orders, we will perform several integrations 
by parts in the integrals appearing in the WKB expansion of the stress tensor, and use the fact that in 
dimensional regularization the integral of a total derivative vanishes [2l|. We will sketch here the calculations, 
leaving the details to the Appendix A. 

The zeroth and second adiabatic orders of the expectation value of the stress tensor can be evaluated from 
Eqs. ([201), (EH) and ([25]). They have been computed in Ref.[H[, and are given by [H]: 

<T W ) (0) = [ (26a) 

(Tu) (V> = — / dxx = — w fc (26b 

4 (2tt)' 1 1 J n - 1 

(r "'" )(2) = 4(2^)»-i / ^^7-{i-(/-- + 2) 2 +^("-l)(/-n + 2)+eG w }, (27a) 
V ; 4(27r)»- 1 7o w fc \ 32(n-l) 1 V ' J 32(n - 1) 1 j 



x 



[4/ - f + (n - 2) 2 ] + ^-[4/ - 2/ 2 + nf + (n - 2)(n - 4)] - - n + 2) + £G u j . 



Performing several integrations by parts and discarding surface terms, one obtains [13j : 



,n — 7i 



IT \(°) — T /oo\ 

(V) 4 ( 27r )n-l 7 °' 

w"'- e >w{ t (,-t^ «44 <29) 

where we have denoted by (i=0,l,2) the integrals given in Table |TJ Eqs. (12"5| and (12T)]) show explicitly 
that the zeroth and second adiabatic orders can be absorbed into a redefinition of the bare cosmological and 
Newton constants respectively. 

The fourth adiabatic order can be computed by following the same procedure. Starting from Eqs. (f^O]) 
and (pH) for (T vv ) and (Tn), we use the adiabatic expansions given in Eqs. (f^Sj) and (j25[) to arrive at the 
following expressions for the fourth adiabatic order of these expectation values: 

r) n™~ n 

(T w ) < 4 > = 4C " (2 ^ (ra _ 1) [aiW'W 2 + a 2 W"W + a 3 W 4 + a 4 H' 2 ] , (30a) 

(Tll>(4) = 4(7(2^"-!) [/?lH ^ 2 + + ^ + /34?< ' 2 + m "' ] - (3 ° b) 

The coefficients and (3i are given in terms of integrals that involve Qk and its derivatives. As anticipated, 
one can find relations between them in order to show explicitly the geometric structure of the adiabatic 
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h) = J Q dxx i co k 


(« + 3) „ o 

l - r°° dx x 2 " fc 


(n-3) 


(»+5) , 2 
7 6 Jo dX 5)5 d3 ,4 


(n+1) o i 




(n-3) 


f" + 7 ' .5-2 

J 8 -J dx—^-^ 




{n + 7) 3 2 2-2 

r _ r°° 2 d d 


/4 = /rdx ^_|_£^ 


*9 Jo ax 7 ^ da;2 < 



TABLE I: Explicit expressions for 1^. To obtain these integrals we have made the change of variables x = k 2 /C and 
we have defined u k = ujk/VC. 



stress tensor. To carry out this procedure, it is convenient to express all these coefficients in terms of the 
integrals I{ defined in Table [IJ For example, let us consider the coefficient ct2, 



f f + °° dxx^ 



Ct-2 = , 

64 Jo 



[4(n- !)(£-£„) + /] 



Taking into account the definition of / (Eq. ([24|) ) wc have 

-./' 



+ °° dxx^ 



G4 



dxx~ 



-f 2 - (31) 



L 



+°° dx. 



("-3) 



2 

3 7. 
2 



- ^3 — o ^3, 



and therefore 



dx 

j 2 _ — -|- _ 

dx dx 3 ,/n dx 
i(n-4)(n-6)/ 3 + |/4. 



(n-1) 2 ^ iy r , (n-4)(n-l) r 7 4 

Q2 = c H H • 

4 V 6 / 1440 160 



+ / 3 



(32a) 



(32b) 



(33) 



Using a similar procedure, outlined in Appendix A, it can be shown that the coefficients are related in 
such a way that the regularized expression for {T^ u )^ takes the form 



(34) 
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where H^J and HjfJ are the tensors given in Eq. (fT4")l . and 



r (3) 



4(2tt)("- 1 ) 
f2„_!/i fi -" f J 3 (n-6) 



(»~4) 
360(?7.-l) 

/ 4 



4(2tt)("- 1 ) 1,1440(71-3) (n-4)(n-3) 
3 (32 J 5 (n - 5) + 60 J 6 - 210 J 7 ) ) 
+ 1440 (n - l)(n - 2)(rc - 3)(n - 4) J ' 



160(n- 

(» + 2) 
4(77-2) 



9(77 + 26) 
1440(77 - 1) 



(35a) 
(35b) 



where all the integrals are shown in Table Q] From Eq. (|3"4"|) we see that the fourth adiabatic order of the 
energy momentum tensor can be absorbed into redefinitions of the bare gravitational constants a, /3, and 7 
of the SEE m 



V. RENORMALIZATION OF THE STRESS TENSOR IN 1+1 DIMENSIONS 

Knowing the dependence with k of the 2j— adiabatic orders, one can use power counting arguments to 
see which of the adiabatic orders contain divergences. In 1 + 1 dimensions, it is simple to check that for 
any of the dispersion relations we are considering, even for the usual one, the zeroth adiabatic order is 
divergent, while the higher orders are finite. This suggests that only a redefinition of the cosmological 
constant would be required to absorb the divergences and, then, one would define the renormalized stress 
tensor as (T A1 ^) re „ = (!),„) — (!),„) However, as will be explained below, this naive argument is incorrect. 
For the usual dispersion relation, it is known that to obtain the renormalized stress tensor it is necessary to 
subtract not only the zeroth adiabatic order, but also the second one As for the usual dispersion relation 
I2 = 0, from Eq. (|2"9"l) we obtain 

O tl n—n 1 

2tt l? 6' 

The point is that the Einstein tensor (which results from the variation of the would be Gauss-Bonnet 
topological invariant at n=2) vanishes as 77 — > 2: G^ v oc n — 2 (see Eq. (fT5"]) ). Therefore, from Eq. (f3"6")) we 
see that the second adiabatic order of the 77-dimensional {T^ u ) is finite in the limit 77 —> 2. However, when 
written in terms of G^ v an explicit pole at 77 = 2 appears. As the pole must be absorbed into the bare 
constant before taking the limit 77 — > 2, a redefinition of the bare Newton constant is also required. In fact, 
the contribution of the second adiabatic order yields the well known trace anomaly in the case of a massless, 
conformally coupled field (£ = 0) [f|: 

(Tp ren = -(T[:)^ = JL. (37) 

On the other hand, for a generalized dispersion relation such that uj^ ~ k 2r with r > 1, the integral 
Ii is finite, but the pole appears explicitly in Eq. ([29|) multiplying I2 (which is also finite). Note that the 
subtraction of (T^)^ 2 ' can be consistently done provided that it is proportional to and hence the pole 
can be absorbed into the bare Newton constant. Since the renormalization prescription must be equivalent 
to a redefinition of the bare constants in the effective Lagrangian of the theory, we conclude that the second 
adiabatic order should also be subtracted in this case. Therefore, the renormalized stress tensor is given by 

Q 



2 +1 " £ +°<"- 2 > 



irp \ irp \ _ irp \ (0) _ irp v(2) 



(38) 
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As a consistency check of this renormalization prescription, let us consider a massless field with £ = £2 = 
in de Sitter space-time, C(rj) = a 2 /r] 2 , and a particular dispersion relation of the form lo\ = k 2 +2b\\k / C '(17) , 
with fen > 0. Wc will compute the trace anomaly by taking the limit in which the dispersion relation tends 
to the usual one (&n — > 0). If we only subtract the zeroth adiabatic order, we have 

Pfl - TO™ = -i jf a (1 - (^Irf - f ) ( 39 ) 

The modes of the field satisfy (see Eq.([7])) 

With the substitution s = (2bu) 1 / i a~ 1 / 2 krj and introducing the constant A = a(26 11 ) -1 / 2 , the equation 
becomes 

^? + (A + S 2 )x fe = 0. (41) 

The particular solution of this equation that satisfies the normalization condition JHJ) and tends to the 
adiabatic mode of positive frequency for \s\ — > 00 (77 — > —00) is given by [l^ : 

X fc ( S ) = ^ r y D_ (1 ^ )[(1 -i)-], (42) 

where I? is the parabolic function @, s = krj/y^X and A = aj \J2b\\. After changing variables and some 
algebra we get 

R f°° , a f „, x VA 

where /(A, s) = fc|xfc(s)| 2 and i? = 12a -2 . A numerical evaluation of this integral gives, in the limit b\\ — > 0, 

(TP (TZ)<® - -A (44) 

As for the case of the usual dispersion relation (fen = 0), the trace of the stress tensor has an anomaly. 
However, the numerical value does not coincide with the usual one (it differs by a sign). Therefore, if we 
subtract only the zeroth adiabatic order, there is a discontinuity in the renormalized stress tensor as &n — > 0. 
This discontinuity disappears if we also subtract the second adiabatic order. Indeed, from Eq.(|29p we find, 
near n = 2, 

<?P (2) = ((2 - n)h + I 2 ) M 2 - (45) 

As I\ is finite for non vanishing coefficient bn, the first term does not contribute to the trace in n = 2. On 
the other hand, I2 is independent of b\\, and an explicit evaluation gives 

< T ;> (2) =-l|- («) 



(Tit) - (T£>c°> = — j o dss* i /(a, S ) - 2 y x ; = [> , 



So, combining Eqs. (|38|) . (|44[) and (|46|) we see that the usual trace anomaly (|37| is recovered in the limit 
611 ^0, 

(r;), e „ - (r;) - <t;>(°) - (r;)< 2 ) - A. (47 ) 

For the sake of completeness, we compute the trace of the stress tensor (T), renormalized according to 
the prescription in Eq. ([551) . for all values of bu (we recall that the limit m — > has to be taken at the end 
of the calculations). In Fig. [T]we have plotted the trace T as a function of A = a/\f2bi\. In this figure we 
see that as A increases (bn decreases) the trace approaches its anomalous value. 
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FIG. 1: The renormalized trace of the stress tensor T (with the prescription in Eg. (|38[) ) normalized to its anomalous 
value as a function of A = a/V26n. These numerical results correspond to the massless scalar field with £ = and 
propagating in a two dimensional de Sitter background. 



VI. RENORMALIZATION OF THE STRESS TENSOR IN 3+1 DIMENSIONS 



In this Section, we will show that in 3 + 1 dimensions the fourth adiabatic order can be consistently 
subtracted, and that this must be done for any of the dispersion relations of the form given in Eq.([B]). 

In 3 + 1 dimensions, by power counting one can show that for u> 2 ~ k 2r , with r > 4, all contributions 
of second or higher adiabatic orders are finite. The divergences come only from the zeroth order terms 
contained in (T^). In the cases io\ ~ fc 6 and ui 2 ~ k A 1 though no fourth order divergences appear, second 
order terms include divergent contributions, which suggests that no terms quadratic in the curvature would 
be necessary in the SEE. However, the situation in 3 + 1 dimensions is similar to that in 1 + 1 23] ■ By using 
the definition of Hj^J given in Eq.(|12p. we can rewrite Eq. (|3"4")) as 



Bo 



(n-4) 



(48) 



As the coefficient B3 does not vanish in four dimensions, an explicit pole at n = 4 appears, which must be 
absorbed into the bare constants of the effective gravitational action ^j. 

If the fourth adiabatic order is not subtracted, there would be a discontinuity in the limit in which the 
dispersion relation tends to the usual one. This means that the renormalized stress tensor would contain 
non vanishing trans-Planckian contributions, even when Mq — *■ 00. We illustrate this point by computing 
the trace of the energy momentum tensor for a massless field with conformal coupling £ = 1/6 and uj 2 = 
k 2 + 2b 'n/c 4 / 'C '(77) in de Sitter space-time, C{rf) = a 2 /r] 2 (as we have done for the 1 + 1 dimensional case). 
From Table [J one can see that for the usual dispersion relation the only divergent integral appearing in the 
fourth order is ^3 (see Eqs. (|34p and (|3"5)) ). So, near n = 4, we have 



1 



2880tt 2 

which is a well known result @, [l^. Therefore, the usual trace anomaly is given by 

R 2 I 

V ^' ren K 34560tt 2 240^ 2 a 4 ' 

On the other hand, for b\\ > all integrals are finite and can be explicitly computed. We find 



(TP 



M\(4) 



1 



1207r 2 a 4: 



(49) 



(50) 



(51) 
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which is —2 times the usual trace anomaly. Note that the result is independent of b\±. 




FIG. 2: The renormalized trace of the stress tensor normalized to its anomalous value as a function of A = aj \J2b\\ . 
The numerical result corresponds to the massless scalar field with £ = 1/6 and uj 2 . — k 2 + 2&nfc 4 /C{rj) propagating 
in a four dimensional de Sitter background. 



The trace of the unrenormalizcd stress tensor with the subtraction of the zeroth and second adiabatic 
orders can be written in the form 

(TP <2£><°> <2£>« = [ + °°d SS s \f(X,s) 



2ir 2 a 4 J Q [ JK ' ' 2(A + s 2 )V2 8(A + s 2 ) 5 / 2 16(A + s 2 ) 7 / 2 

By means of a numerical evaluation, as bn — > 0, we obtain 



(52) 



{Tp _ (T , )( o) _ (T , )(2 ) ^ (53) 

Therefore, the usual trace anomaly is recovered in the limit &n — » only when the fourth adiabatic order is 
also subtracted: 

(Tftren = (T») (T^ (T^ (T^ (54) 

The behaviour of the renormalized trace T as a function of A = a/y/2bn is shown in FigfJl where we see 
that it approaches to its anomalous value as the dispersion relation tends to the usual one. 



VII. SELFCONSISTENT DE SITTER SOLUTIONS 



As an application of the results presented in the previous sections, here we will study the SEE in de 
Sitter space-time, C(rj) = a 2 /rj 2 , taking into account quantum effects of free scalar fields with a generalized 
dispersion relation (i.e., including the backreaction of the quantum fields on the spacetime dynamics). 

It is well known that in de Sitter space-time the expectation value of the energy momentum tensor {T^ v ) of 
a free scalar field with the usual dispersion relation is proportional to the metric <7 M „ in the de Sitter invariant 
states [26j . As they are related by Bogoliubov transformations with momentum independent coefficients, 
these states form a one complex parameter family. Here, we will point out that for a generalized dispersion 
relation there also exists a one complex parameter family of states for which (T^) cx g llUl which is just 
the tensorial structure of the cosmological constant term. We will then analyze the relation between the 
de Sitter curvature (R — 12a -2 ) and the cosmological constant (A) for a dispersion relation of the form 
lo 2 = k 2 + 2b 11 kyC{r 1 ). 
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In de Sitter space-time, the field modes \k satisfy (see Eq.©) 



drj 2 y k\u ' ??2 

where (l 2 = m 2 + n(n — 1)(£ — £„)/a 2 . With the substitution s = krj/yX, this equation can be recast as 

— + Iw (s) + — — J x fe = 0. (56) 

where u> 2 (s) = Loj,(r))/k 2 is a function of fc and r\ only through the variable s. Let /(s) and g(s) be two 
independent solutions of this equation. Then, a field mode can be conveniently expressed as 

Xk (s) = ^f(s) + ^g(s). (57) 

We now choose the coefficients = A and = B to be momentum independent. Thus, defining ip(s) = 
Vkxk(s), the normalization condition {8j becomes 

^{s)^-(8)-^(s)r(s)=iVx. (58) 

By introducing these particular solutions (Xfc( s ) = ii>{s)/Vk) into Eqs. (fl7|) and (|T8|) . and rescaling the 
integration variable one can show that p = {T vri )/C and p = (Tn)/C are time independent. Therefore, for 
the corresponding states we have (T^) ex g^, provided that {T^ v ) is covariantly conserved. It is remarkable 
that this one-parameter family of states exists for any dispersion relation. 

If we choose the particular state of the family that reproduces the WKB solution as \s\ — > oo, the 
divergences in the stress tensor can be absorbed into the bare gravitational constants, and the SEE can be 
recast as 



G>w + Ag^v = »irNG N {Tp,) rm , (59) 

where we N is the number of scalar fields and we assumed that the renomalized values of a, (3 and 7 vanish. 
Using that in de Sitter space-time i? M „ = and (T /J ^) ren = Tg^ v /A (where T = (Tjf) ren ), we have 

+ A = 2ttNG n T. (60) 

4 

Therefore, the de Sitter metric is a consistent solution of the SEE even when the backreaction of scalar 
quantum fields with a generalized dispersion relation is included. 

For a massless field with £ = 1/6 and uj' 2 . = k 2 + 2bnk 4 /C{rf) {jl 2 = 0), Eq. (f5"6")) can be solved exactly. The 
solution for the modes are given in Eq.|@U) (where s = {2b 11 ) 1 / A a - 1 / 2 k'q = kr]/VX). Note that as b n -> 
these field modes tend to the standard Bunch-Davies modes, defining a generalized Bunch-Davies vacuum. 
To get the renormalized trace T of the stress tensor we subtract up to the fourth adiabatic order. An integral 
expression for the unrenormalized trace with the subtraction of the zeroth and second adiabatic order is given 
in Eq. ([52)) . while in Eq. tjSTj) we have the fourth adiabatic order. Note that the function / = T/R 2 depends 
on only one free parameter buR = 6/ A 2 . Therefore, we have 

A(R) = — + 2nNG N T = — + 2nNG N R 2 f(b 11 R). (61) 

We evaluate f(buR) numerically for different values of b\\R. The relation in Eq. (f6"Tj) is shown in Fig. [3J 
where we have also plotted the classical relation A = R/4. The results are presented for &h?ti 2 j = 1 and 
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FIG. 3: The relation in Eg. (|61fl between the curvature of the de Sitter space-time R and the cosmological constant 
A for two values of &n: In solid line &n = m~ 2 and in dashed line bn = 10 4 m~ l 2 (where m v i is the Planck mass), 
with N = 1 (on the left) and N = 100 (on the right). In each case, the classical relation A = i?/4 is shown by a 
dotted line. These results correspond to N massless scalar fields with £ = 1/6 and u\ = k 2 + 2buk 4 /C(rj) in four 
dimensions. Note the different scales in both axes. 

10 (where m p i is the Planck mass) for N — 1 and N — 100 (for intermediate values of the parameters the 
results lie in between). From these figures, we can see that for any of the values of the parameters the general 
features are common. There is no selfconsistent solution for large values of the cosmological constant A. The 
value of A above which there is no more selfconsistent solution becomes larger as bn decreases. When A is 
small there are two selfconsistent solutions: one is near the classical solution A = R/4 while the other has a 
positive curvature R larger than 4A, even for negative values of A, which in general is far from being in the 
semiclassical regime. These results are similar to those obtained in Ref.[27| for the usual dispersion relation. 

VIII. AN EFFECTIVE INITIAL STATE? 

It has been suggested in the literature that the trans-Planckian effects could be taken into account in a 
low energy effective field theory with usual dispersion relations, by considering a generic "initial" state for 
the modes of the quantum field when they leave the sub-Planckian regime [HI . There is some debate about 
whether the trans-Planckian corrections to the power spectrum of primordial fluctuations can be consistently 
reproduced from a suitable choice of the initial state [H, [l?], [HI • Moreover, the choice of the initial time for 
imposing the initial condition is a nontrivial problem. In fact, as was pointed out in Ref . [30| . an inadequate 
choice could lead to artificial oscillations in the power spectrum. On the other hand, within this description 
it is difficult to quantify the backreaction effects [l2j . 

We will discuss now whether the backreaction effects on the spacctime metric can be taken into account by 
considering an arbitrary initial state in the low energy effective theory. If we adopt the usual renormalization 
prescription in the effective theory, the divergences in (S , |T At „|S') can be absorbed into counterterms of the 
gravitational effective action if the state \S) coincides with the adiabatic state up to the fourth order (see for 
instance [H, Hil)- Hence, the (3k coefficient of the Bogoliubov transformation that relate the mode function 
V>f (corresponding to the state \S}) to the Bunch-Davies one ipj? D , 

^ = a k ^ D +M^ D \ (62) 

is required vanish faster than k~ 2 as k — > oo. In previous works (T5 , IT?! |29| it has been shown that the de 
Sitter invariant states belong to a one parameter family, which is related to the Bunch-Davies vacuum by 
constant Bogoliubov coefficients. Therefore, the only renormalizable state which is de Sitter invariant is the 
Bunch-Davies vacuum. 
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In Section IVHi we have pointed out that for any generalized dispersion relation, a family of states for 
which (T^„) oc g^ v also exists (as for the de Sitter invariant states in the usual theory). Thus, by choosing 
the member of the family that tends to the adiabatic mode of positive frequency for |s| — > oo, one can obtain 
a generalized Bunch-Davies state for which the renormalized stress tensor is proportional to the metric 
tensor. Therefore, there are states for which the expectation value of the stress tensor of a scalar field with 
a generalized dispersion relation cannot be reproduced by any choice of the initial quantum state of a field 
with the standard dispersion relation, unless other renormalization scheme is employed [l8| . In other words, 
while can obtain selfconsistent de Sitter solutions for any generalized dispersion relation, in the standard 
theory this is possible only for a single quantum state. 



IX. FINAL REMARKS 



In this paper we have presented a complete analysis of the renormalization procedure for the semiclassical 
Einstein equations, in theories in which the quantum scalar fields satisfy modified dispersion relations. This 
work generalizes the adiabatic renormalization for quantum field theory in curved spaces developed in the 
seventies to theories containing higher spatial derivatives of the matter fields. 

We have shown that, even though power counting suggests that for this class of theories it would be 
enough to renormalize the cosmological and Newton's constants, in 3 + 1 dimensions a consistent procedure 
also involves the subtraction of the fourth adiabatic order. Therefore, it is also necessary to include terms 
quadratic in the curvature into the gravitational action in order to absorb the divergences of the expectation 
value of the energy momentum tensor. This subtle point was missed in our previous work fl3j , and clarified 
in the 1 + 1 dimensional case in Ref . [l4| • 

We obtained regularized expressions for the stress tensor up to the fourth adiabatic order, and showed 
explicitly the geometric nature of its divergences. We also computed the trace of the renormalized tensor 
in de Sitter space for the 'would be conformal field' in the standard theory, and recovered the usual trace 
anomaly in the limit M c — > oo. 

We have also shown that, whatever the dispersion relation, there exist a family of quantum states for 
which the mean value of the stress tensor is proportional to the metric. One member of this family (the 
'generalized Bunch-Davies' vacuum), is renormalizable, and therefore it allows the existence of selfconsistent 
de Sitter solutions, as explicitly shown for massless, conformally coupled fields. These solutions are not 
present for arbitrary states in the standard theory (they only exist for the usual Bunch-Davies vacuum), and 
therefore this particular trans-Planckian effect can not be simulated by modifying the quantum state in an 
effective theory with the usual dispersion relation. 



Appendix A: Regularization of (T M „) 



(4) 



In this Appendix we outline the procedure by which we derived Eq. (|34[) . starting from the unrenormalized 
fourth adiabatic order of the stress tensor. The idea is to find relations between the different integrals that 
appear in the WKB expansion of the stress tensor, in order to reveal its geometric nature. The relations can 
be found by performing integrations by parts and by using that in dimensional regularization the integral of 
a total derivative vanishes [2l|. 

^From Eqs. (f2T))) and (j2"Tj) . with the use of the WKB expansion, the fourth adiabatic order of (T vv ) and 
(Tn) can be recast as 



<r w > (4) 



Q„_l£t f 



4C(2tt)(* 



{2) wln 2 



~,3 



-/< + (n -!)(£- £n) 3/- 



[ {2) w 2 ] 2 - ^w 2 'n 

(n-2) 



t +2(n -!)(£-&) 



(63) 
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<T n ) (4) 



n—n z'+oo 



4C(2tt)("- 1 ) J 8w£(n-l) 



+ (2) w fc 2 'w 

+ (2 % fe 2 H 2 



{ (2) wf' [/ + 4(n - - £„)] ~ [ (2) ^ 2 ] 2 (2« - 5 - 3/) 



2(n-l)(f-e„)(6/-(n-l))-^ + |/ 



fn — 91 

6K-^)(n-l)/+V/+/ 2 



(£ - £„)(« - 1) -(n - l)(n - 2) + 9/ 2 - 6/ - 3(n - 1)/ + + 



(n-7) 



(64) 



where the function / = /(.t) is defined in Eg . (f2~4"|) . Here we have used Eg. ([23)1 to write in terms of 

and its derivatives. By using in addition that 



(2) wf = ^[4« - £«)(n - 1) + /] + ^[/ 2 - 6/ + 4(n - 2)(n - - &)] + ^[2/ - //], (65a) 
(2) wf = ^[/ 2 - 6./ + 4(n - 2)(n - - £„)] + ^[/ 2 - 8/ + 4(n - 2)(n - - &)] (65b) 

'1/4 T_/2a_y/ 7-/'" 

+ ^[/ 2 + // - 2/'] + —[12/ - 5//] + — [4(£ - £ n )(n - 1) + /], 

we write (T w )( 4 ) and (Tn)( 4 ) in the form given in Eq. (|3"0"|) . where the coefficients on and Pi can be expressed 
as a linear combination of integrals of the form 



HH' 



rnnls 



(66) 



with m,n,l,s integer numbers. For example, the coefficient <xi considered in the text (Eq.(|3ip) is given by 

1 



a2 = -\{n - 1) 2 (C - Cnfh m - UJ- 



1000 — 77^2000) 

64 



where I3 = Joooo is defined in Table [H and from Eq. (|32p we identify 

(n - 4) T 



•^1000 

"^2000 



_(n-4)(n-6)J 3 + |-f4. 



It is straightforward to show that 



(67) 

(68a) 
(68b) 



a 4 = ~2 a 2- 

To find a relation between a.\ and ai requires a little more work. The coefficient ol\ is given by 

1 



In order to express it in terms of the integrals 1^ in Table [J it is useful to observe that 



J: 



2000 



X 2 



/ 



1 



— —^iooo + 2 
(n-4) 



x duil. 
Cj\ dx 



x 2 o(a; fc 
cD 4 . (ix 



JlOOO + g^OlOO, 



(69) 



*i = -r(n - 1) 2 (£ - f») Jiooo - r^(n - 1)(€ - £n)[3J 2 ooo - 2J ioo] - -^[3J 3 ooo - 4J 1100 ]. (70) 



(71) 
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where the last equality follows after performing an integration by parts and discarding the surface term. 
Thus we have 

3J2000 — 2 J0100 = (n — 4) J10007 (72) 

which is one of the combination of integrals appearing in Eq. (|70p . In a similar way, the integral J3000 can 
be recast as 



'3000 



'2000 



J)2 dx 

x 2 du) k j, 2 
<D| dx 



-J2000 H g ^2000 + g^llOO: 



from which we have 



3J3000 - 4Jnoo = (n - 4) J 200 o- (74) 
Inserting Eqs. (JBSaJ), {72) and {TjJ into Eq.([7D)). we find 

on = ~(n - 1) 2 (C - Cn) 2 (n - 4)/ 3 - ^(n - - £„)(n - 4)J 1000 - ^(n - 4)J 20 oo, (75) 
and comparing this coefficient with a2 (given in Eq. (|67|0 we see that 

ai = ^y^a 2 . (76) 

Therefore, with the use of Eqs. p4a[) . (|14c[) . (f6"9") and (fTB") . the component {T m )^ can be written as 

r(!) qow(3) 



/T \(4) _ ^"-1 M" " J a 2H nn i2H m 
\ W ' : 4 (2^)n-i | ( n -l) 2 ~ ( n -l)(n-2)(n-3) 



(n- 10)(n-2) 
" 3 32 " 2 



(77) 



Since from Eqs. (f3"Tj) and (f6"8"| the coefficient a2 can be expressed in terms of the integrals Ii, Eq. (|35a[) for 
_Bi follows straightforwardly. 

To show that the same equation (|77[) is satisfied by the component (Tn)^, with i/^ replaced by if^ , 
we start by noting that 

ft = -T-^-rraa, (78) 
(n- 1) 



(e-6 



\2 



A = * ^ (n-1) 



-- + (n-4) 



-^3 + ~^~~7p~ [6J2000 — 4Jqioo — ^1000] (79) 



6^/3000 — 8J1100 — ^2000] 1 



64(n-l) 

& = ^ ~ 6l)2 (n - 1) [9Jiooo + h] + ($ ?/ n) t 9j2 "°" ~ 6 Joi °" + Jl0 "°] ( 8 °) 
o lb 
1 

[9 J3000 — 12Jnoo + ^2000] • 



128(n- 1) 
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With the use of Eqs. (jtJgaj) . (|?2"1) and (75) we have 



/h = 

04 = 



-- + (n-4) 

1 3(ti-4) - 
4 4 



As, 
/9s. 



Thus (Tu)( i can be expressed in the form 



(Tn) 



(4) 



4(27r)™^ 1 I (n-1) 2 

03 



(|(„-4)(n-6)-|) 



C 



^ + i V 1 (" 2 -13n + 28))/3 5 



7T 



(81) 
(82) 

(83) 



By the same procedure of integrating by parts and discarding surface terms, and after a lot of algebra, we 
find 



0s = ^0i 



(n-5) 



2(n-2)(n-3) 



B 3 , 



(84) 



where the coefficient S3 is given in Eq. (|35a[) . Finally, we arrive at Eq.([34]) by inserting these last relations 
into Eqs. G3 and ([83]>. 
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